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We investigate the transmission of electrons between conducting nanoribbon leads oriented at 
multiples of 60° with respect to one another, connected either directly or through graphene polygons. 
A mode-matching analysis suggests that the transmission at low-energies is sensitive to the precise 
way in which the ribbons are joined. Most strikingly, we find that armchair leads forming 120° 
angles can support either a large transmission or a highly suppressed transmission, depending on 
the specific geometry. Tight-binding calculations demonstrate the effects in detail, and are also used 
to study transmission at higher energies as well as for zigzag ribbon leads. 
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I. INTRODUCTION 

Graphene, a two-dimensional honeycomb lattice of car- 
bon, is one of the most interesting new low-dimensional 
materials to have become available in the laboratory in 
the last few yearsi. When undoped, the low energy 
physics of this system is dominated by two Dirac points^ . 
The wavefunctions associated with states near them are 
described by spinors, whose amplitudes represent the 
probability density to find electrons on either of the two 
honeycomb sublattices. Because the fermion spectrum is 
gapless, these spinors have well-defined helicity, leading 
to an absence of backscattering from impurities--"^. The 
observed metallic behavior of undoped graphene is likely 
to be a manifestation of this suppressed backscatteringi. 

Among the interesting and potentially useful proper- 
ties of graphene is the prospect of "tailoring" its elec- 
tronic properties by cutting it into ribbons of well-defined 
widths along various symmetry directions^i^i^ii. Re- 
cent experimental work*'^'"'^^ has confirmed the possibil- 
ity of tuning transport gaps of graphene ribbons via their 
widths, although the quality of these ribbons is not yet 
high enough to be usefully compared with the expected^ 
width dependence of ideal ribbons. In applications one 
generically needs to join such ribbons together as inter- 
connecting wires or elements of a device. Understanding 
the transport through such junctions is then important in 
designing graphene geometries with desirable behaviors. 
This the subject of our study. 

Changing the direction of electron currents would pre- 
sumably be an important aspect of any graphene-based 
circuit. In both quantum wires and electromagnetic 
waveguides it is known^^ that the transmission through 
bends depends on the detailed nature of the bend geome- 
try. This raises the prospect that the conductive behavior 
of a graphene junction may differ substantially from the 
properties of its individual nanoribbon leads. 

In the simplest situation, the nanoribbon leads meeting 
at a junction are of the same type (armchair or zigzag), 
restricting the bend angles to either 60° or 120°. The 



latter case is particularly interesting due to the behavior 
of the low energy states near the Dirac points under 60° 
rotation. These states may be constructed, within the 
k • p approximation, from products of the exact wave- 
functions at the Dirac points, which vary rapidly in real 
space, with slowly varying envelope functions^i. The ro- 
tation induces a transformation on the fast component of 
the wavefunction which exchanges both the valleys and 
the sublattices. As a result, states near the Dirac points 
are nearly orthogonal to their 60 degree rotation in a 
confined geometry. Since a ribbon with a 120° junction 
may be viewed as 60° deflection of the electron trajectory, 
one might expect a generically suppressed transmission 
through 120° bends. 

Our studies show that, while this naive reasoning is 
sometimes borne out, in general the transmission through 
such bends is not universal and depends critically on the 
details of the junction. In this paper, we focus on ge- 
ometries in which 120° bends in armchair nanoribbons 
are realized either by a "kink" , or by attachment to tri- 
angular or hexagonal central regions. We focus mostly 
on the energy region \Ef\ < Ee, E^ being the band 
edge of the first excited band of the nanoribbon leads, 
where there is only one channel available for conduction. 
In this low energy region, we know that the eigenstates 
of graphene nanoribbons may be understood within a 
continuum approximation (the Dirac equation) when ap- 
propriate boundary conditions are adopted^, so that a 
description involving matching of these wavefunctions at 
jmictions becomes possible. As we will discuss in detail, 
these geometries produce quite different transmission be- 
havior at low energies. This suggests the prospect of 
using different types of junctions to tailor transmission 
through a set of graphene nanoribbons. 

To understand why the geometry is crucial, at low 
energies we may adopt wavefunctions for armchair rib- 
bons obtained from solutions to the Dirac equation^, and 
consider how they might be appropriately matched at 
a junction. We focus on armchair nanoribbons whose 
widths are chosen so that they are metallic^. A "mode- 



matching" procedure may be formally developedi^ to 
compute conductance properties of the system. This be- 
comes particularly simple when only the lowest energy 
transverse modes are retained - the "single-mode approx- 
imation" (SMA). 

In the simplest case, the 120° junction (see Fig. [T]), 
we shall see that the in the zero energy limit, incoming 




FIG. 1: Schematic diagram of a 120° junction, illustrating the 
coordinate systems for the two ribbons [{xi,yi) and (3:2,2/2)]. 
The surface at which the two nanoribbons are joined is the 
solid vertical line. 

and outgoing modes can be matched up perfectly at the 
junction, so that at low energies the electrons are maxi- 
mally transmitted. By contrast, passage through a short 
length of zigzag ribbon (Fig. [21) involves intermediate 

Zigzag 




FIG. 2; Schematic diagram of a 120° armchair junction, in 
which electrons pass through a length of zigzag ribbon. Coor- 
dinates for the three sections illustrated. Wavefunctions must 
be matched along the two solid lines connecting the armchair 
ribbons to the zigzag ribbon. 

transverse states that are strongly localized to the edges, 
which become orthogonal to those in the armchair leads 
very close to zero energy. This implies blocked transmis- 
sion in a narrow range of energies near zero. Transmission 
through an equilateral triangle with two attached arm- 
chair leads may be viewed as a special case of this class 
of geometries. 

A closely related geometry is transmission through an 
equilateral triangle with three leads. Here the system 
may be constructed from three appropriately cut arm- 
chair leads as illustrated in Fig. \3\ As we will show 
in detail for this case, rapid oscillation of the fast com- 
ponent of the wavefunction makes the transmission very 
sensitive to the precise way in which these armchair leads 
are connected to the triangle, so that one may obtain ei- 
ther large or vanishingly small conductance in the low 



energy limit. This sensitivity to the geometry is ubiqui- 
tous for such graphene nanostructures, suggesting that a 
wide variety of conductance properties can in principle 
be engineered into very similar geometries. 




FIG. 3: Schematic diagram of an equilateral triangle with 
three armchair leads, in which the triangle is taken as an 
"endcap" for the vertical lead. Coordinates for the three sec- 
tions illustrated. Wavefunctions must be matched along the 
two solid lines of the triangle. 



The SMA thus leads one to expect qualitatively dif- 
ferent conductances for 120° armchair junctions at low 
energies, depending on precisely how they are joined. 
We have tested these expectations using a tight-binding 
model^^ for the ribbons and the regions joining them. 
Again, the simplest case is the 120° junction between two 
armchair nanoribbons. An example of a specific connect- 
ing geometry and its associated low-energy conductance, 
computed in a tight-binding model for ribbons of con- 
ducting width, is illustrated in Fig. 31 In agreement 
with the SMA result, one obtains almost perfect trans- 
mission at low energies, nearly to the bottom of the first 
excited transverse subband energy. 

Fig. m illustrates a typical result for two leads joined 
through an equilateral triangle. As suggested by the 
SMA, the conductance is now suppressed near zero en- 
ergy, but only over a very narrow range. Tight-binding 
studies of analogous geometries in which armchair leads 
are connected by short zigzag ribbons give similar results 
at low energies. One may also compute the conductance 
through 120° when a third lead is added to the equilateral 
triangle. The results of such calculations are illustrated 
in Figs, inland [71 In the former case, for which the leads 
are attached to the triangle at their widest cross-sectional 
widths, one finds a highly suppressed conductance at low 
energy. The latter case has the leads attached to the tri- 
angle at their narrow cross-section, and the resulting con- 
ductance is large at low energy. We shall see both these 
results may be understood from the SMA, as special cases 
of the geometry illustrated in Fig. I3l Thus, the mode- 
matching procedure appears to offer a useful framework 





FIG. 4: Conductance per spin of 120-degree bends in arm- 
chair nanoribbons having N transverse channels. Geometry 
for N^8 illustrated in (d). 



FIG. 6: Conductance per spin of triangles with 3 armchair 
nanoribbon leads having N transverse channels, with ribbons 
connecting to the triangle along their wider cross-section. Ge- 
ometry for N=8 illustrated in (d). 



for understanding conductance in such geometries. 






(d) Geometry for N 



FIG. 5: Conductance per spin of triangles with 2 armchair 
nanoribbons leads having N transverse channels. Geometry 
for iV=8 illustrated in (d). 




FIG. 7: Conductance per spin of triangles with 3 armchair 
nanoribbon leads having N transverse channels, with ribbons 
connecting to the triangle along their narrower cross-section. 
Geometry for A''=8 illustrated in (d). 



The tight-binding calculations also allow us to examine 
the conductance at higher energies, as well as to consider 
other geometries. Here we summarize a few results from 
such studies, and give further details later on. (1) Trans- 
mission studies through hexagons with armchair leads re- 
veal low-energy conductances that are also suppressed at 



low energies for 120° transmission. However, in this case 
the suppression is quite dramatic throughout the lowest 
subband. Conductances per spin through 60° and 180° 
is a finite fraction of e^/ft. in the lowest subband. (2) At 
higher energies, we find in all of the geometries consid- 
ered that the transmission as a function of energy approx- 



imately follows the density of states for bulk graphene, 
provided the ribbon widths are large enough. However, 
very close to the energy of the van Hove singularity we 
find a strong suppression or enhancement of the conduc- 
tance depending on the angle between the leads. (3) 120° 
zigzag ribbon junctions have a more complicated evolu- 
tion with increasing ribbon width than their armchair 
cousins. Most notably, very close to zero energy the con- 
ductance oscillates between large and small values as the 
ribbon width is incremented by a single unit, reminiscent 
of recent results for p — n junctions of zigzag ribbons^^. 
The remainder of this article is organized as follows. 
In Section H, we discuss the mode matching analysis 
in more detail, and explain how the SMA leads to the 
expectations described above. Section III is devoted to 
describing the numerical methods used to compute the 
conductance through the various geometries in the tight- 
binding model. In Section IV we provide further details 
of our numerical results. Finally, we conclude in Section 
V with a summary of our results. 



II. MODE-MATCHING ANALYSIS 

The conduction properties of electron systems in which 
current is injected into and removed from a region with 



a defined shape and potential, through infinite leads 
with known cross-section, can be understood by exploit- 
ing their analogy with electromagnetic waveguides^^. A 
conceptually simple approach to their analysis is to di- 
vide the system into components where the wavefunc- 
tions may be computed with appropriate boundary con- 
ditions, and then "stitch" the wavefunctions together at 
the boundaries by matching their amplitudes and, in 
the case of wavefunctions controlled by the Schroedinger 
equation, their derivatives. Such calculations are made 
analytically tractable by employing the single-mode ap- 
proximation (SMA), in which only the lowest subband 
of the external leads is retained, which produces qualita- 
tively and often quantitatively good results provided one 
works away from scattering resonances of the systemic. 
A similar strategy may be employed to understand our 
low energy numerical results for systems with armchair 
leads. 

The simplest case to consider is the 120° junction be- 
tween two armchair ribbons illustrated in Fig. [TJ The 
geometry is divided into regions 1 and 2 with correspond- 
ing wavefunctions V'*-^^ and "0^^^ and coordinate systems 
{xi,yi), {x2,y2)- At low energies one may write down 
approximate forms for the ribbon eigenstates^. For mo- 
mentum p these have the form 
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where the upper (lower) entry represents the amplitude 
on the A (B) sublattice, W is the ribbon width, the 
Vi — (xi,yi) represent the positions of lattice points on 
the two ribbons, K = (-|j,0), K' = (|f,0), and a is 
the lattice constant. The value oipx must be chosen such 
that the total amplitude at the edges of the ribbons van- 
ishes, hence Px -^ Pn comes in quantized values^. For 
metallic ribbons, the lowest subband satisfies Pn=o = 0. 
These wavefunctions have energy e — wf|p|, where vf 
is the speed of electrons near the Dirac points. Due to 
particle- hole symmetry, we restrict our analysis to e > 0, 
setting e = Ep for the determination of the conductance 
at zero temperature. 

A general wavefunction in which current is injected 
only from the left in the lowest subband (0 < £ < vppi) 
may be written in the form 

v^(i) = Ao<u^-f5:s(^)<j,_,„ (2) 

n 

0(^) = j:B'Mlpy (3) 

n 

In this expansion, it is implicitly understood that for 
modes where vppn > e, one replaces ipy with ±Kj, in 



Eq. [T] such that the wavefunctions appearing in Eqs. [2] 
and [3] are evanescent as one moves away from the junc- 
tion. The S^'^ coefficients are determined in terms of Aq 
by matching the wavefunctions, for both sublattices, on 
the solid line shown in Fig. [TJ The conductance per spin 

is then given by G = ^(iB^'^lVl^oP)- 

To carry out this procedure, one must first specify a set 
of matching conditions at the joining surface that guar- 
antees continuity of the wavefunctions and the current 
across the junction. One possible choice is illustrated in 
Fig. [HJ where it is now convenient to change notation 
slightly and refer to wavefunctions tp, tp' and coordinates 
(x, y), {x' , y') in regions 1 and 2, respectively. Equating ij) 
to "0' on the joining line (open and closed circles) guaran- 
tees continuity of the wavefunctions. Matching currents 
across these junction points can be more complicated be- 
cause this in general involves products of wavefunctions 
on either side of a bond. However we can greatly sim- 
plify the latter matching condition by anticipating the 
SMA, for which we will use the wavefunctions of Eq. [U 
which vanish on the open circles in Fig. [5] for the low- 
est subband. Thus one only need match the currents on 
the bonds connecting the closed circles to the triangles. 




5' 4' 3' 2' r 



FIG. 8: Diagram detailing the joining surface between leads 
for the simple 120° bend, as well as the labelling scheme for 
X and x' coordinates. Matching the wavefunction on the dark 
circles and triangles accounts for current continuity across the 
joining surface. 



This may be accomplished straightforwardly by matching 
the wavefunctions on the triangles as well. We focus on 
the zero energy transmission and take the Py —^ limit. 
The resulting matching conditions may now be written 
explicitly in the form 



(4) 



45(a; = 3,y) = 47(x' = 3',2/') 

4^J{x = 6,y) = 4f{x' = 6',y') 
4^^{x^7,y)^4''J'ix'^7',y'), 
i;^,l\x^8,y)^4^^'{x' = 7',y'), 
4^J{x = 9,y)^4^J'{x' = 9',y') 



where the x and x' labels are defined in Fig. [H] Simple 



phase factors may be added to these matching conditions 
to generalize them for the case Pj, ^ 0. Note that in 
writing these equations we have equated A sites of the left 
lead with B sites of the right. This may be understood 
if one constructs the junction by starting with a single 
ribbon, excises an equilateral triangular region from the 
center, and joins the two resulting ribbons to form the 
junction. In doing this one finds that at the junction, 
A sites are indeed matched up to B sites. Note that no 
explicit interchange of the K and K' valleys is needed 
(as is the case in 60° rotations) because the real space 
coordinate system has also been rotated. 



Eqs. [3] are a realization of setting Eq. [5] to Eq. [3] on 
the joining surface. To proceed, we wish to represent the 
wavefunctions on the matching points in an expansion 
in terms of wavefunctions of the form in Eq. [TJ For- 
mally, this is accomplished by multiplying these equa- 
tions by 4'w n' (2^7 y); &^^d then integrating (x,y) on the 
joining surface, i.e., summing over the points where the 
wavefunctions have been matched. (Note p' is chosen 
such that Vp(p[^ ~^ P'y) ~ ^^0 This results in a set of 
equations relating the ^0 and i?*^*^ amplitudes. A sec- 
ond set of equations may be generated by multiplying 
the matching equation by V'l' „' i^: v) and integrating on 
the joining surface. The two sets form an in principle in- 
finite dimensional matrix equation that relates the i?^*^ 
coefficients to Aq. 



Carrying out this procedure is vastly simplified by 
adopting the SMA, in which only the lowest transverse 
mode, Pn = 0, is retained in the matrix equationi^. To 
demonstrate the perfect transmission at low energy in the 
junction illustrated in Fig. [Tl it is convenient to consider 
the reflection amplitude in the SMA. This is proportional 
to the integral 



Mo,o(p,) ^ d\ 4^1;^(a;(A),z/(A))42^(a;(A),y(A)), 



(5) 
where A parameterizes the joining surface. Note that 
in the limit py — > 0, there is no actual y (jj') depen- 
dence in ^Q^{x,y) (i/igQ (a;',y')). From Fig. [5] one may 
see that the positions denoted as x{') = n(') demarcate 
increments of length a/2. The meaning of the formal 
expression (Eq. [5]), using Eq. [U then takes the form 
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That -Mo.o(Py) vanishes in the limit Py ^> Q indicates 
an absence of backscattering as e ^ 0, and hence perfect 
transmission in this hmit. One may also compute the 
corresponding overlap on the joining surface for trans- 
mission and confirm that it has a magnitude of unity. 
Deviations from this are of order {pyWY , so that these 
become significant when py ~ l/VF, which occurs at an 
energy of the same order as the bottom of the first excited 
subband. Thus for energies well below this, we expect the 
transmission to be very close to unity. This behavior is 
confirmed by the tight-binding calculations. 

This behavior seems dramatically different than the 
naive expectation discussed in the Introduction, that the 
interchange of the K and K' valleys might lead one to ex- 
pect a 60° deflection of the electron trajectory to be sup- 
pressed. The mode matching procedure however demon- 
strates that only the overlap on the joining surface need 
be considered, and because this involves a small sub- 
set of lattice points, destructive interference between the 
rapidly oscillating parts of the wavefunction may not be 
realized. 



J 



A second example of this procedure may be consid- 
ered for the geometry illustrated in Fig. [51 in which 
two armchair leads are joined at the two solid lines to 
a short length of zigzag nanoribbon. Traversal through 
the equilateral triangle with two leads may be thought 
of as a special case of this geometry, with the short- 
est possible zigzag ribbon. Approximate wavefunc- 
tions *&p„,p„(a;3, ys) for the zigzag ribbon region may be 
developed^. These are more complicated than the arm- 
chair forms, in that both p„ and Py vary as e varies, so 
that the transverse wavefunctions vary with e even within 
a single subband. At energies close to zero, this varia- 
tion becomes quite pronounced in that the wavefunctions 
become highly localized at the zigzag ribbon surfacesii. 

One may develop an explicit expression for the trans- 
mission amplitude for this geometry, within the single 
mode approximation, in terms of the overlap integrals on 
the two junctionsi^. The result is proportional to the 
product of the overlap integrals on each of the joining 
surfaces, Ngl{py,p'y), with 



NgoiPy^Py) = 



dAi<^;(a;(Ai),y(Ai))$(Jp,(x(Ai),y(Ai)) 
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In these integrals, Ai and A2 parameterize the left and 
right surfaces in Fig. [21 zLpo{e) is the transverse mo- 
mentum from which the lowest zigzag transverse state is 



made up^ii^, and e^ 
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The 



portant observation in this case is that, at low energy, 
the states of the zigzag ribbon $ become confined to the 
surfaces, with a length scale ^(e) which vanishes rapidly 
at low energy in the continuum description^. Since the 
lowest transverse state of the armchair ribbon remains 
spread throughout the ribbon cross-section, one may see 



(«) 



-^00 ^ ^/^(e)/^^! and the resulting conductance will be- 
have as G ^ [^(e)/^^]^. This means the conductance per 
spin is is suppressed near zero energy, but can rise to a 
value of order e^/h once £ is above the range of energies 



where zigzag ribbons support surface states. The result- 
ing conductance is suppressed in a narrow range near 
zero energy. 

As a final example, we consider the three lead equilat- 
eral triangle geometry. The SMA has a form identical 
to the case of the zigzag ribbon described above, with 
the overlap integrals now performed on surfaces joining 
armchair ribbons. Unlike the above two examples, the 
surface is oriented at different angles with respect to the 
cross-sections of the two ribbons, as illustrated in Fig. 
[3l In this case the two integrals Nq q whose product 
is proportional to the transmission amplitude involve a 
product of wavefunctions whose fast components vary at 
different rates as one moves along the joining surface. 



The resulting overlaps are very sensitive to the precise 
way in which the leads are joined to the triangle. 



FIG. 9: Diagram detailing the joining surface between the 
lower lead and the upper left lead, for an equilateral triangle 
with three leads, and the leads attached to the triangle at the 
wider point in their cross-section. 




Fig. [9] details a joining surface for the geometry of 
Fig. [6l Note in this case that the ribbons are joined 
along zigzag edges, so that the matching of both wave- 
functions and currents is easily accomplished. Denoting 
the wavefunction in the lower lead by iPqI {x, y) and the 

one in the upper left lead by V'o^p i^' , y'), with fi ~ A,B, 
we find the matching conditions for py — to have the 
form 



4^0^ {x = 1, y) = 4fix' = 1', y') , Vfo^ (x = 2, y) = ^'^'{x' = 2', y') 
4^ {x = 3, y) = 4f(x' = 2', y') , ^^^^ {x = 4,y) = ^'^'{x' = 3', y') 
4-^ {x = 5, y) = 4f{x' = 3', y') , 4"^^ {x = 6,y)^ 4f{x' = 4', y') 



(11) 
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Again, for py = there is no actual y (y') dependence 
in 4oQ{x,y) (V'ao (^''2/'))- The overlap between the in- 



coming state from the bottom lead and the outgoing state 
from the upper left lead, using Eq. [U has the form 
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The upper line in Eq. [T2]is due to the overlap of the wave- 
functions on the A sites, while the lower line comes from 
the B sites. Multiplying out the square brackets gener- 
ates terms which are either independent of the integer n 
or oscillate in n with period 3. One finds that the non- 
oscillating terms from the upper and lower lines precisely 
cancel. The remaining rapidly oscillating terms vanish in 
the sum provided the maximum value of n is a multiple 
of 3, and in any case give vanishing contribution as the 



ribbons become wide {W — > oo). The cancellation of the 
non-oscillating terms indicates a complete destructive in- 
terference between incoming and outgoing waves for the 
two arms of the triangle at low energy, as one might have 
naively supposed. With no overlap at the joining surface, 
the conductance should vanish at zero energy. 

The other simple geometry for joining the ribbons to 
the triangle is detailed in Fig. [101 In this case the match- 



ing conditions take the form 
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and the corresponding overlap sum is now 



-'^0,0 






exp 



-E 



W 



27r 

-i — n 

3 



27r 
exp [ I — n 



2tt 



I Att , lA f Att, 1, 

exp ( *y('^+ 2^ ^°^P i^*Y'-"+ 2^ 



47r 
exp I i—n 



exp 



—2 cos ( — ) + 2 + ( oscillating terms) 



An 

-i — n 

3 



(14) 



r 



We see the slight shift in positions of the B sites where 
the wavefunctions are matched in Eq. [T3] leads to an 
extra phase factor, such that the A and B overlaps no 
longer cancel. Thus Nq q is relatively large in this case, 
and we expect a correspondingly large transmission. 

We next turn to our numerical studies, which, as dis- 
cussed in the Introduction, essentially confirm the expec- 
tations of the SMA. 



III. MODEL AND NUMERICAL METHOD 

Our calculations arc based on a simple tight-binding 
model of graphcne, in which only nearest neighbor hop- 
ping is included. Formally the Hamiltonian may be writ- 
ten as 



H = -t 



E 



(|R.)(R,| + |R,)(R,|). (15) 



Here t is the nearest-neighbor hopping matrix element 
and {RiRj} denotes bonds on a honeycomb lattice with 
boundaries defined by the polygon of interest and/or the 
attached leads. For the purposes of this section only we 
use the energy unit t — 1. An example of a hexago- 
nal scattering region with attached leads is illustrated in 
Fig. 111! We consider only armchair and zigzag nanorib- 
bons for our leads, although other periodic edges may 
be considered within our method. We wish to calculate 



conductances through various pairs of leads in this geom- 
etry. Below we discuss the procedures used to evaluate 
the conductance per spin G from Green's functions of the 
Hamiltonian iJ, under the assumption of time-reversal 
symmetry. 

Conceptually, we divide the lattice into three regions: 
a central region C, a "left" lead L and a "right" lead R, as 
illustrated in Fig. [T^lfor the case of armchair nanoribbon 
leads. We compute the conductance between the leads 
L and _R, and any remaining leads in the geometry are 
considered part of the central region C . 

Conductance at zero temperature is computed as a lin- 
ear response using the Kubo formula^", which relates con- 
ductivity to a current-current correlation function. For 
the two-lead conduction problem, the relevant current 
operators correspond physically to the total current in 
the region L flowing in the incoming direction and the 
total current in the region R flowing in the outgoing di- 
rection. We denote these operators as Jl and J_r, whose 
precise deflnition we elucidate below. Charge conserva- 
tion implies that the current flux down a given lead in 
a time-independent, steady state is the same regardless 
of where in the lead it is "measured" . For our purposes, 
it is most convenient to "measure" the current precisely 
at the two boundaries where the the respective leads join 
the central region. These boundaries are each illustrated 
in Fig. [12] as a pair of dashed lines. 

These considerations motivate the following defini- 



boundary 





boundary 



C 







FIG. 12: Conceptual division of the geometry into regions C, 
L, and R, illustrating the lattice point labeling scheme at the 
boundaries, which are denoted by pairs of dashed lines. This 
example is a symmetric armchair ribbon with A*' = 4, and the 
boundary is chosen such that Q = 5. 



FIG. 10: Diagram detailing the joining surface between the 
lower lead and the upper left lead, for an equilateral triangle 
with three leads, and the leads attached to the triangle at the 
narrower point in their cross-section. 




FIG. 11: Hexagon with 6 leads 



tions: we label the lattice points immediately to either 
side of this boundary as r^ ^ '' for lattice points belonging 
to region L{R) and s^ for those belonging to the cen- 
tral region. Here /z runs over integers 1 to Q, where Q is 
the number of bonds traversing the boundary. Typically 
Q is the same as the number of channels N, but Q > A^ is 
possible in certain configurations. This labeling scheme 
is illustrated in Fig. 1121 The current operators may be 
explicitly defined as 



'L{R) 



K~^)T.i\''^''^){^''"' 



p=l 



,LW^^r^(i..) 



(16) 
The Kubo formula for conductivity leads to the well- 



known22 transmission formula for the conductance per 
spin 



G = ^MtLRtLR. + i^jRL) 



(17) 



where iLfl(_RL) is the transmission amplitude for states at 
the Fermi energy Ep from the region L to R {R to L) and 



tLR. 



"RL- 



In our context, the Kubo formula involves 



the operators J^ and Jn, and the resulting transmission 
amplitudes in pT)) are Q x Q matrices. 

To evaluate the transmission amplitude, the relevant 
retarded Green's functions are Q x Q matrices whose 
components are given by 






K\{Ef 



IT] 

irj 






4) 



4) 



(18) 
(19) 
(20) 
(21) 



Here _ff^, H^, H'-^ represent the restriction of the Hamil- 
tonian to regions L, i?, and C, respectively (with no hop- 
ping across the boundaries), and a and (3 can be L or R. 
Ideally, 77 — > 0^, but for numerical calculations it is nec- 
essary to choose a small 77 > which in effect becomes the 
energy resolution of the computation. Our calculations 
are carried out with r] = 10~^i. 

With these definitions, we have, in abbreviated nota- 
tion for the Q X Q matrices, 



'4^g[R.L]vL 



tLR. — 



1/2 
L 



(22) 
(23) 



Schematically, this formula shows that the transmission 
amplitude is given by the propagator from the left to 
right side of the central region, with the vl{r) veloc- 
ity matrices^^ normalizing the nanoribbon states to unit 
flux. The problem is now reduced to computing Q'^^^i 
and Q . Standard gluing formulas^^ for non- interacting 
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Green's functions can be used to obtain Q[R,L] from 
these two. 

For our geometries, Q^^^^ is the Green's function of 
a semi-infinite nanoribbon at its termination. One may 
directly evaluate the Green's function of a single nanorib- 
bon unit cell by inverting the matrix (Ep+irj—h), where 
h is its lattice Hamiltonian. We then rapidly extend from 
the unit cell to a nanoribbon segment of length 2' through 
I successive length-doubling steps via the aforementioned 
gluing formulas. We find that the Green's function on one 
termination of the long nanoribbon segment becomes an 
accurate approximation for the semi-infinite ribbon when 
2' ~ t/f]. This can be verified by substituting the numer- 
ical result into a Dyson equation satisfied by the exact 
Green's function. Such a Dyson equation may be easily 
derived for any such semi-infinite periodic structure. 

For Q"-^ we first suppose that the geometry has only two 
leads and thus the region C is finite. One approach is to 
perform the matrix inversion of {Ep + it] — H^). How- 
ever, significant computational savings are possible when 
the Green's function is only required at the boundary-^. 
In the case of more than two leads, G'^ can be found by 
gluing the Green's functions for the "passive" leads in the 
calculation to the Green's function of the finite scatter- 
ing region. This procedure correctly accounts for current 
which is drained away by the extra leads. 

We first tested our numerical techniques on the 
straightforward case of infinite nanoribbons, whose band- 
structure is well- under stoodl. Nanoribbons may be char- 
acterized by the number N of conducting channels. The 
integer N gives two related properties of the ribbon: 1) 
the minimum number of severed bonds required to break 
the ribbon into two disconnected pieces (as in Fig. [T^]) 
and 2) the maximum possible value of G, the conduc- 
tance per spin, in units of the conductance quantum 
e^//i. Armchair nanoribbons have associated with them 
two lines of fictive lattice points (just outside the ac- 
tual ribbon edges) separated by a width La (see Fig. 
1131) on which the wavefunction vanishes. When the rib- 
bons possess a reflection symmetry through the center, 
N = L — 1, and otherwise N = L — 1/2. In either case, 
the ribbon is semiconducting (i.e. there is gap in the 
spectrum around zero energy) unless sin(47rL/3) = 0. 
Thus, symmetric armchair nanoribbons are metallic only 
for the series A^ = 2, 5, 8, . . . and the asymmetric ribbons 
for N = 1,4,7,10,.... Zigzag ribbons are metallic for 
any value of N. 

We have compared the conductance of these various 
types of ribbons as a function of Ep with their band- 
structure. We find, as expected, a contribution to G{Ep) 
of e'^/h for each band present at Ep, with the excep- 
tion of the flat bands2i in symmetric armchair ribbons 
at energies ±i. We also verify the linear dispersion of 
the metallic band in armchair nanoribbons with veloc- 
ity vp — \/3ta/2, as well as the maximum G values and 
metallic N sequences discussed above. 






FIG. 13: Symmetric armchair, asymmetric armchair, and 
zigzag ribbons with A*' = 4. Dashed line cuts through four 
bonds of each ribbon, the minimal number required to sever 
each into two disconnected pieces. Heavy vertical lines indi- 
cate rows of fictive lattice points along which the wavefunction 
vanishes for the armchair nanoribbons. 



IV. NUMERICAL RESULTS 

A. Armchair Lead Systems: 120° Junctions and 
Triangles 

Our low-energy result for the simple junction for arm- 
chair leads at low energy was discussed in the Introduc- 
tion, and is illustrated in Fig. U) One clearly sees that 
the tranmission is essentially perfect, with a small sup- 
pression just below the first excited subband energy. The 
behavior appears to be well-explained by the SMA. The 
120° junction also naturally occurs in an equilateral tri- 
angle with two leads, as illustrated in Fig. |5(d)[ As has 
already been explained, at low energies one finds suppres- 
sion of the conductance in a narrow range around zero 
energy. 

We discuss more fully the example of the three lead 
triangle geometry, illustrated in Fig. [6{d). Because of 
the symmetry of this geometry, the two point conduc- 
tance is the same for any pair of the three leads. Fig. 
[Til illustrates the conductance for three different system 
sizes as a function of Fermi energy Ep over the entire 
bandwidth. One may see that the overall shape of the 
conductance curve roughly tracks the density of states 
for bulk graphene, with peaks at the van Hove singu- 
larities Ep = t, where t is the hopping matrix element. 
The bumps and wiggles around this are due to changes 
in the number of conducting channels as the Fermi en- 
ergy is increased, as well as to quantum interference in 
the scattering regionii. 

The low energy region of conductance for different sys- 
tem sizes may also be considered. In Fig. [15] we blow 
up the low-energy region for three different system sizes. 
The cases A^ = 13 and A^ = 15 as expected reveal no 
conduction at low energies, since there are no conducting 
states to carry current through the leads. The metallic 
state for A^ = 14, by contrast, allows a finite conductance, 
but one may see its actual value is remarkably small at 
very low energy. One may examine this behavior with 
increasing A", and not surprisingly the pattern of near 
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{h)N=lO 




{c)N=15 

FIG. 14: Conductance per spin through an equilateral triangle 
with three armchair leads of sizes (a) N=5, (b) N=10, and 
(c) N=15. 



suppression for system sizes of the form N — 3M + 2, 
and complete suppression for other sizes, repeats itself. 
When viewed as a function of Ep/Ee, which fixes the 
position of the opening of the first excited subband as 
N becomes large, the conductance in the lowest subband 
tends to a roughly parabolic shape, very small but re- 
maining finite away from Ep = as N —f oo. This 
represents the continuum limit, in which the width W of 
the ribbons remain finite, and the lattice constant a is 
taken to zero. The increase from zero of the conductance 
as the energy increases from zero is in agreement with the 
results of the SMA, but we note that the vanishing con- 
ductance for these widths requires the joining geometry 
illustrated in Fig. [9l As discussed in the Introduction, a 
joining geometry of the form illustrated in Fig. 1101 leads 
to a non- vanishing conductance (Fig. [7]), in agreement 
with the SMA analysis. 



B. Transmission Through a Hexagon 

We next consider the case of a hexagon with six at- 
tached leads, as illustrated schematically in Fig. [TT] The 
conductance per spin is illustrated for the full bandwidth 
in Figs, [m [T71 and [TS] Here we must specify the angle 
between the two leads upon which the measurement is 
made: Fig. [TBI corresponds to a 60° angle between leads, 
Fig. [17] to 120°, and Fig. [18|to 180°. Several remarks 
arc in order. As in the equilateral triangle, the overall 
structure of the conductance follows the density of states 
for bulk graphene. However, as the size of the system 
increases it is apparent that there is a remarkably sharp 
suppression for transmission through 60° and 180°, and 
a strong enhancement for transmission through 120°, at 
the van Hove singularity, Ep/t = 1. The rapid change in 




Ef/E, 
(a)Af=13 



Ef/E, 

(b)7V=14 




Ef/E^ 

(c)Af=15 

FIG. 15: Conductance per spin in the low energy region for 
an equilateral triangle with three armchair leads, (a) A'^=13, 
(b) N^U, (c) iV=15. 



resistance with respect to Fermi energy suggests that this 
phenomenon could in principle be useful as a transistor, 
although the relatively high Fermi energy where it oc- 
curs may require a large electric field to realize. Beyond 
this, it is also noteworthy that the overall scale of trans- 
mission through 120° is significantly larger, and seems to 
grow more quickly with system size, than for the other 
two directions. 




(a)Af=5 




FIG. 16: Conductance per spin for two leads at an angle of 
60°, through a hexagon with six armchair leads, (a) N=5, 
(b) iV=15, (c) iV=20. 

The conductance at low-energy for this geometry is il- 
lustrated in Fig. [THJlor N — 20, for conductance through 
each of the three angles. A suppression of 120° transmis- 
sion is apparent, and in this case is in fact much more 
pronounced throughout the first subband than is the case 
for the equilateral triangle. Indeed, in studies of hexagons 
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(a)7V=5 



{h)N=15 






FIG. 17; Conductance per spin for two leads at an angle of 
120°, through a hexagon with six armchair leads, (a) 7V=5, 
(b) iV=15, (c) iV=20. 



FIG. 19: Conductance per spin for two leads at an angle of 
(a) 60°, (b) 120°, and (c) 180° through a hexagon with six 
A^=20 armchair leads. 





(c)Ar=20 

FIG. 18: Conductance per spin for two leads at an angle of 
180°, through a hexagon with six armchair leads, (a) N=5, 
(b) iV=15, (c) iV=20. 



with only three or two leads, we find the 120° transmis- 
sion to be even more suppressed throughout the lowest 
subband than in the six lead case. This suggests the 
hexagon may be useful in three terminal devices where 
one may wish to employ one lead as a voltage probe with- 
out draining current flowing between other leads. 



C. Zigzag Nanoribbon Junctions 

We conclude this section with a summary of analogous 
results for a 120° zigzag nanoribbon junction. Like arm- 
chair ribbons, zigzag ribbon widths may be characterized 



by the minimum number of broken bonds N required to 
sever it, as illustrated in Fig. [T31 Unlike armchair rib- 
bons, zigzag ribbons are metallic for any widthS.. At low 
energies the current-carrying states have an interesting 
chirality in that left- and right-moving states are asso- 
ciated with different valley indicesiS. It should be em- 
phasized that the association of this discrete index can 
only be made precise in a continuum description. All 
the geometries considered below may always be under- 
stood in terms of semi-infinite ribbons joined together at 
boundaries on the lattice scale. Because the lowest en- 
ergy states of zigzag ribbons are highly confined to the 
edges of the system, a pure continuum description is in- 
adequate even at the lowest energies. Similar physics has 
been noted recently in graphene p — n junctions^. 

This physics is most clearly seen in 120° zigzag junc- 
tions. Fig. [20] illustrates the transmission for this ge- 
ometry for ribbon widths N—6,7, and 8. At very low 
energies, there is a qualitative difference between odd 
and even width ribbons, with the former supporting a 
large conductance at zero energy and the latter a small 
one. Such odd-even behavior also occurs in p — n junc- 
tions, and appears to be related to the fact that edges in 
a zigzag ribbon align when TV is even, but anti-align for 
odd iVMJI. Even at higher energies, but still within the 
lowest subband, the conductance as a function of energy 
appears to change qualitatively from one width to the 
next as N increments by single units. 

The importance of lattice scale physics in this system 
is further made apparent by an examination of the local 
density of states at low energies. This is shown in in Fig. 
[21] for a junction of width N—6 at several energies below 
the first excited subband energy. At low energies, while 
the wavefunction is strongly maximized near the ribbon 
edges, the junction can attain zero or perfect conduc- 
tance, as is the case in (a) and (b), respsectively. As the 
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(d) Geometry for A'^ 



FIG. 20: Transmission through 120° bend for zigzag ribbons, 
(a) iV=6, (b) iV=7, (c) iV=8. (d) Geometry for iV=8. 



wavefunctions throughout the first subband intimately 
involve the lattice scale, any continuum description for 
this system is not likely to be reliable. 



V. CONCLUSION 

In this paper we have studied the conductance of 
various graphene geometries in which the current in a 
nanoribbon is redirected through angles that are multi- 
ples of 60°, which can be accomplished in the honeycomb 
network without introducing lattice defects or changing 
the ribbon type. We focus on the low-energy behavior 
in armchair nanoribbon geometries and find a variety of 
behaviors, including very high transmission for a partic- 
ular realization of a simple 120° junction and suppressed 
transmission for the same angle when there is an inter- 
vening triangle. A mode-matching analysis, within the 
single-mode approximation (SMA), allows one to under- 
stand in a simple way many of these results. With this 
technique we demonstrate that the rapid oscillation of 
the low energy wavefunctions renders the conductance 
through such junctions highly sensitive to the precise 
way in which the ribbons are joined together. Tight- 
binding calculations support the conclusions of the SMA 
at low-energies, and further elucidate the details of the 
conductance at higher energies. 



We also presented numerical results for conductance 
through other geometries. Hexagons in particular showed 
a dramatic suppression of conductance through 120°, and 
further supported a strong enhanced/suppressed conduc- 
tance (depending on the angle between the leads) at the 
van Hove singularity. Zigzag nanoribbon junctions were 




{a.)EF/Ee=0.0l6 



{h)EF/Ee=0.16 




(c)EF/Ee=0A9 



{d)EF/Ee=0.82 



FIG. 21: Local density of states shown in gray scale for 120° 
degree zigzag junction, N — 6. Black denotes zero, and white 
denotes maximal LDOS. The leads modelled in the calculation 
extend away from the junction to infinity. 



also studied, and were found to have a richer behavior 
than their armchair cousins, which is likely to require a 
more microscopic description than is possible in a con- 
tinuum model. 
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